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ABSTRACT 

Here  we  report  on  an  effort  to  include  an  empirically  based 
transition  modeling  capability  in  a  RANS  solver.  Testing  of  well- 
known  empirical  models  from  literature  for  both  attached-  and 
separated- flow  transition  against  cascade  data  revealed  that  the  models 
did  not  provide  enough  fidelity  for  implementation  in  an  airfoil  design 
system.  Consequently,  a  program  was  launched  to  develop  models 
that  would  provide  sufficient  accuracy  for  use  in  an  airfoil  design 
system.  As  a  first  step  in  the  effort,  accurate  modeling  of  free  stream 
turbulence  development  was  identified  as  a  need  for  any  form  of 
transition  modeling  capability.  Additionally,  capturing  the  effects  of 
free  stream  turbulence  on  pre- transitional  boundary  layers  was  found  to 
have  a  significant  effect  on  the  accuracy  of  transition  modeling.  A 
CFD- supplemented  database  of  experimental  cascade  cases  (57  with 
attached- flow  transition  and  47  with  separation  and  turbulent 
reattachment)  was  constructed  to  explore  the  development  of  new 
correlations.  Dimensional  analyses  were  performed  to  guide  the  work 
and  appropriate  non-dimensional  parameters  were  then  extracted  from 
CFD  predictions  of  the  laminar  boundary  layers  existing  on  the  airfoil 
surfaces  prior  to  either  transition  onset  or  incipient  separation.  For 
attached- flowr  transition,  exploration  of  the  database  revealed  a  distinct 
correlation  between  local  levels  of  freestream  turbulence  intensity, 
turbulence  length  scale,  and  momentum- thickness  Reynolds  number  at 
transition  onset.  It  was  found  that  the  correlation  could  be  recast  as  a 
ratio  of  the  boundary- layer  diffusion  time  to  a  time-scale  associated 
with  the  energy-bearing  turbulent  eddies.  In  the  case  of  separated- 
flow  transition,  it  was  found  that  the  length  of  a  separation  bubble 
prior  to  turbulent  re-attachment  was  a  simple  function  of  the  local 
momentum  thickness  at  separation  and  the  overall  surface  length 
traversed  by  a  fluid  element  prior  to  separation.  Both  the  attached- 
and  separated-flow  transition  models  were  implemented  into  the 
design  system  as  point- like  trips. 

INTRODUCTION 

In  axial- flowr  turbomachinery,  the  design  trend  is  toward 
increasing  airfoil  loading  in  an  effort  to  reduce  weight  and  cost  of 
future  systems.  Transition  prediction  is  critical  for  accurate  loss 
predictions  of  high  lift  airfoils,  and  the  full  multi -moded  (Mayle,  1991) 
nature  of  the  transition  process  must  be  considered.  Lakshminarayana 
(1991),  Simoneau  and  Simon  (1993),  Simon  and  Ashpis  (1996),  Dunn 
(2001),  and  Yaras  (2002)  all  provide  detailed  reviews  of  the  state  of 
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the  art  in  predictive  techniques  for  turbomachineiy,  and  they  point  to 
the  need  for  improved  models  for  transition. 

Elevated  levels  of  freestream  turbulence  (7Y/>1.0%)  have  a 
significant  effect  on  pre-transit ional,  or  quasi-laminar  boundary  layers. 
Further  to  that,  it  is  the  authors7  opinion  that  the  quality  of  the  laminar 
boundary  layer  at  transition  onset  must  be  captured  accurately  before 
any  form  of  transition  modeling  can  be  employed.  Therefore,  it  is 
important  to  capture  accurately  the  field- wise  development  of 
freestream  turbulence  quantities.  To  that  end,  the  ability  of  the  k-m 
turbulence  model  of  Wilcox  (1998)  to  predict  the  development  of  Tu 
was  validated  against  the  experimental  data  of  Ames  (1994). 
Additionally,  an  accurate  technique  for  modeling  the  effects  of  Tu  on 
laminar  boundary  layers  within  the  frame w'ork  of  the  k-6)  model  was 
developed. 

In  testing  against  cascade  data  it  was  found  that  open- literature 
models  for  attached  and  separated-flow  transition  were  not  sufficiently 
accurate  for  implementation  in  a  design  system.  Consequently,  an 
effort  was  launched  to  develop  new  correlations  for  attached-  and 
separated-flow  transition.  A  dimensional  analysis  was  performed 
considering  all  transition- relevant  quantities  available  within  the 
frame-work  of  a  RANS  simulation  performed  with  a  two-equation 
turbulence  model.  A  database  of  the  resulting  dimensionless  groups 
was  constructed  from  open-literature  and  Pratt  &  Whitney  in-house 
cascade  data.  The  cascade  data  were  supplemented  with  quantities 
based  on  the  aforementioned  modeling  techniques  for  freestream 
turbulence  development  and  its  effects  on  laminar  boundary  layers. 
An  investigation  of  the  resulting  database  enabled  the  development  of 
new  models  for  attached-  and  separated-flow  transition.  The  details  of 
this  process  are  documented  below, 

A  computational -methods  section  will  be  presented  first  with 
discussions  on  boundary'  conditions,  freestream-disturbance 
propagation  and  quasi-laminar  boundary  layers.  Then,  sections 
concerning  attached-  and  separated- flow  transition  are  presented, 
where  reviews  of  the  state-of-the-art  and  current  model  development 
details  are  discussed.  Validation  and  benchmarking  of  the  new  models 
is  presented  in  Part  II  of  this  manuscript. 

COMPUTATIONAL  METHODS 

Steady-state  and  time-resolved  turbine  flowfields  were  predicted 
using  the  3-D,  Reynolds -Averaged  Navier-Stokes  (RANS)  code 
described  both  by  Ni  (1999)  and  Davis  et  al.  (1996).  Numerical 
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closure  for  turbulent  flow  is  obtained  via  the  k-m  turbulence  model  due 
to  Wilcox  (1998).  An  O-H  grid  topology  was  employed  for  all 
simulations,  and  approximately  600,000  grid  points  per  passage  were 
used  for  three-dimensional  simulations  executed  for  this  study 
(without  tip  clearance).  The  viscous-grid  provides  near- surface  values 
ofy+  less  than  1  over  all  no-slip  boundaries  and  gives  approximately  7 
grid  points  per  momentum  thickness  in  airfoil  and  endwall  boundary 
layers.  These  grid  densities  and  spacings  are  consistent  with  Pratt  & 
Whitney  standard  work  criteria  defined  by  Praisner  and  Rangwalla 
(200!)  to  capture  thermal  fields,  surface  heat  transfer  and  transition- 
related  streamwise  gradients  in  gas  turbines*  The  code  is  accurate  to 
second-order  in  space  and  time  and  multi -grid  techniques  are  used  to 
obtain  rapid  convergence.  Constant  temperature,  constant  heat  flux 
and  adiabatic- wall  thermal  boundary  conditions  are  available  and  were 
employed  when  appropriate*  The  inlet  boundary  conditions  used  for 
each  simulation  are  described  below. 


The  Development  of  Freestream  Disturbances 

Freestream  disturbances  have  been  shown  to  play  an  important 
role  in  transition  (See  e*g.  Abu  Ghannam  and  Shaw  1980,  and  Mayle, 
1991)*  In  order  to  develop  a  RANS-based  transition-modeling  system, 
some  form  of  freestream  turbulence  modeling  must  be  implemented  in 
the  solver.  Two-equation  turbulence  models,  which  provide  modeling 
of  the  convection,  dissipation,  and  diffusion  of  freestream 
disturbances,  are  currently  the  state  of  the  art  in  turbomachinery  design 
systems.  However,  the  ability  of  any  two-equation  model  to  capture 
accurately  the  development  of  freestream  turbulence  quantities  in 
prototypical  turbomachineiy  flowfields  should  be  demonstrated  before 
empirical  modeling  based  on  freestream  disturbances  can  be  applied. 

Two-equation  turbulence  models  present  the  user  with  two 
quantities  that  must  be  specified  at  the  inlet  of  a  simulation*  1)  The 
turbulence  kinetic  energy,  k,  is  directly  proportional  to  the  local 
turbulence  intensity  and  velocity  and,  when  simulating  experimental 
data,  is  typically  set  based  on  measured  levels*  2)  The  dissipation 
parameter  should  be  set  to  match  the  experimental  decay  rate  of 
freestream  turbulence  when  it  is  derivable  from  published  data*  When 
it  is  necessary  to  estimate  inlet  boundary'  conditions  for  two-equation 
turbulence  models,  there  are  a  number  of  techniques  that  may  be 
employed. 

When  simplifications  for  zero  pressure-gradient  steady  flow  are 
applied  to  the  k-m  equations  (Wilcox,  1988)  one  can  obtain  the 
follow  ing  relations  for  the  freestream  development  of  k  and  m: 


*(*)=*,„ 
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where  x  is  streamwise  distance,  kin  and  mtn  are  inlet  quantities,  and  f/„ 
is  the  freestream  velocity.  If  the  measured  turbulence  decay  rate 
upstream  of  the  cascade  is  known,  equations  1  and  2  can  be  used  to 
solve  for  the  inlet  values  of  k  and  m *  The  predicted  decay  rate  of 
freestream  turbulence  obtained  from  equations  1  and  2  varies  with  x  to 
the  -0*62  power,  and  this  falls  within  the  range  of  -0.60  to  -0.68 
reported  by  Baines  and  Peterson  (1951)  and  Hinze  (1975)* 

If  the  freestream  turbulence  decay  rate  upstream  of  the  test 
section  is  not  know  n  then  there  are  three  other  means  of  deriving  the 
inlet  value  of  the  specific  dissipation  rate,  mm.  Firstly,  if  the  decay  rate 


is  not  reported  for  a  configuration  with  grid  generated  turbulence,  and 
either  the  grid  location,  or  the  bar  size  of  the  grid  is  known,  then  the 
decay  rate  can  be  estimated  quite  accurately  by  using  the  following 
relation  which  is  similar  to  one  from  Baines  and  Peterson  (1951): 

/  \ -0.65 
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Here  d  is  the  cross-stream  dimension  of  the  bar  elements  that  comprise 
the  grid*  If  the  experimentally  estimated  dissipation  rate  is  available, 
the  inlet  value  for  m  can  be  estimated  based  on  the  relation  from 
Wilcox  (1998): 
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where  £is  the  measured  dissipation  rate  and  C^=0.09. 

Finally,  as  a  third,  less-preferrfcd  methodology,  the  authors  have 
found  that  the  following  relation  from  Wilcox  (1998)  gives  a 
reasonable  estimate  for  m  at  the  inlet  in  cases  where  the  length  scale 
for  an  experimental  configuration  is  reported  as  well  as  the  inlet 
turbulence  level: 


Q)-- 


(5) 


where  X  is  the  measured  length  scale  and  u ’  is  deduced  from  the 
freestream  turbulence  level.  Use  of  equation  5  for  grid-generated 
turbulence  typically  results  in  inlet  values  for  m  close  to  those  obtained 
with  equation  1. 

Using  these  techniques  for  deriving  k  and  m  boundary'  conditions, 
comparisons  were  made  between  CFD  simulations  with  the  k-m  model 
and  the  measured  freestream  development  of  turbulence  quantities  in 
the  vane  cascade  from  Ames  (1994)  and  Ames  and  Plesniak  (1995), 
In  his  reports,  Ames  provides  detailed  measurements  of  turbulence 
quantities  within  the  passage  of  the  C3X  vane*  In  one  case  freestream 
turbulence  was  generated  with  a  passive  grid,  providing  a  turbulence 
intensity  of  8%  nominally  at  the  inlet  to  the  cascade*  This  test  case  is  a 
good  compromise  between  configuration  accuracy  and  completeness 
of  measured  flow'  quantities  for  this  type  of  benchmarking* 
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Figure  1 .  Measured  and  predicted  distributions  of  freestream 
turbulence  around  the  C3X  airfoil.  Data  is  from  Ames  (1994). 

Figure  1  is  a  comparison  of  the  measured  and  predicted 
turbulence  levels  around  the  airfoil.  Also  shown  is  the  distribution  of 
Tu  obtained  if  '‘frozen’*  turbulence  is  assumed  within  the  passage. 


2 


Frozen  turbulence  assumes  that  only  the  airfoil  velocity  distribution 
affects  the  local  Tu  level.  This  method  results  in  over-predictions  of 
local  levels  of  Tu  of  up  to  40%  on  the  pressure  side  and  in  the 
stagnation  region  while  giving  good  estimates  on  the  suction  side. 
Overall,  better  agreement  is  obtained  with  the  freestream  modeling 
supplied  by  the  k-m  model.  Although,  experimentally,  the  isotropy  of 
the  freestream  turbulence  is  significantly  reduced  by  the  potential  field 
of  the  airfoil,  the  evolution  of  the  simulated  isotropic  turbulence  is 
within  ±0.5%  of  the  measured  levels  based  on  the  stream- wise 
component  of  turbulence  on  both  the  pressure  and  suction  sides  of  the 
airfoil 

For  accurate  multi -row  modeling  of  multi-stage  turbomachinery 
turbulence  fields,  it  is  also  important  to  capture  accurately  the 
development  of  the  turbulence  length  scale/dissipation  rate.  Figure  2 
is  a  comparison  of  the  predicted  and  experimentally  estimated 
dissipation  rates  around  the  airfoil  tested  by  Ames  (1994).  While  the 
agreement  between  CFD  and  data  is  not  as  good  here  as  in  Figure  1,  it 
is  shown  in  Part  II  of  the  present  work  that  this  level  of  fidelity  is 
indeed  adequate  for  making  multi -row  predictions  of  attached-  and 
separated-f!owr  transition.  Further,  predicting  the  development  of 
freestream  turbulence  through  a  turbine  passage  as  in  Figure  1  is  itself 
dependant  upon  accurately  calculating  the  rate  of  turbulence 
dissipation  as  the  flow  proceeds  downstream.  Again,  further 
validation  of  the  ability  of  the  k-m  model  to  predict  the  development  of 
freestream  disturbances  in  multi-stage  environments  will  be  presented 
in  Part  II  of  this  paper. 
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Figure  2.  Comparison  of  measured  and  predicted  turbulence 
dissipation  around  the  C3X  airfoil  Data  is  from  Ames  (1994). 


Modeling  Pre-Transitional  Boundary  Layers 

Moss  and  Oldfield  (1992)  concluded  from  their  experimental 
study  of  turbulence  level  and  length  scale  on  heat  transfer  in  laminar 
boundary  layers  that,  “accurate  prediction  of  heat  transfer 
enhancement  due  to  freestream  turbulence  is  not  possible  using  the 
turbulence  level  alone,”  Considering  Reynolds’  analogy,  it  seems  that 
accurate  predictions  of  quasi- laminar  thermal  and  momentum 
boundary  layers  can  not  be  obtained  without  attention  to  turbulence 
length  scale  (i.e.  specific  dissipation)  as  well  as  intensity. 

The  importance  of  capturing  turbulence  intensity  and  length-scale 
effects  on  laminar  boundary  layers  was  emphasized  in  the  recent  report 
of  Boyle  ef  ah  (2003)  who  developed  a  model  that  considers  both 


turbulence  level  and  length  scale  to  increase  turbulent  viscosity  above 
zero  in  laminar  regions  to  account  foj  what  they  refer  to  as  “turbulence 
enhancement/"  Boyle  and  coworkers  found  a  “noticeable” 
improvement  in  spatially  averaged  laminar-region  heat  transfer 
predictions  using  their  model,  but  no  model  tested  produced  accurate 
local  heat  load  levels.  Similarly,  Roach  and  Brierley  (2000)  point  to 
the  importance  of  modeling  the  effects  of  turbulence  level  and  length 
scale  on  pre-transitional  boundary  layers.  However,  the  authors 
assume  that  the  integral  quantities  of  quasi-lam  mar  boundary  layers 
are  the  same  as  the  equivalent  purely  laminar  (7«=0)  boundary  layer. 
It  will  be  demonstrated  that  this  assumption  is  not  consistent  with 
experimental  data.  Sharma  et  al.  (1988)  reviewed  experimental 
evidence  that  turbulence  from  upstream  rows  imparts  a  significant 
influence  on  the  pre-transitional  boundary  layers  on  turbine  airfoils. 
Additionally,  the  importance  of  capturing  the  effects  of  freestream 
turbulence  on  laminar  boundary  layers  has  been  identified  in 
experimental  studies  of  laminar  heat  transfer  rates  by  Ames  (1995)  and 
Van  Fossen  et  al.  (1994). 

The  phenomenon  of  a  quasi-laminar  boundary  layer  is 
demonstrated  in  the  cascade  heat  transfer  data  of  Arts  et  al.  (1990), 
Figure  3  is  a  plot  of  measured  heat  transfer  distributions  from  Arts  et 
al.  ( 1 990)  with  nominal  inlet  freestream  turbulence  levels  of  1%  and 
6%.  Additional  flow  conditions  for  both  data  sets  in  Figure  3  were 
Rec=lxlG6  and  M^-Q.77,  where  Rec  is  Reynolds  number  based  on 
true  chord.  Heat  transfer  augmentations  in  the  leading-edge  and 
pressure-side  regions  of  the  airfoil  of  up  to  approximately  50%  are 
evident.  Given  that  the  only  difference  between  the  two  experiments 
is  the  turbulence  level,  the  enhanced  energy  transfer  across  the  quasi- 
laminar  boundary  layer  in  the  high  Tu  case  must  be  primarily  due  to 
the  penetration  of  freestream  turbulence  into  the  laminar  boundary 
layer. 
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Figure  3.  Measured  convective  heat  transfer  coefficient 
distributions  from  Arts  et  al.  (1990)  and  CFD  predictions  run 
with  fully  laminar  boundary  layers. 


The  data  from  Arts  et  al.  (1990)  was  employed  to  assess  various 
methods  for  modeling  quasi-laminar  boundary  layers.  Uniform 
surface- temperature  conditions  were  held  for  each  of  the  simulations 
consistently  writh  the  measured  values.  CFD  simulations  (Figure  3) 
were  run  fully  laminar  on  both  the  suction-  and  pressure-sides  of  the 
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airfoil  employing  common  modifications  from  literature  to  simulate 
laminar  regions.  One  such  technique  tested  that  is  commonly 
employed  in  transitional  RANS  simulations  involves  setting  the  eddy- 
viscosity,  juT,  from  the  turbulence  model  to  zero  within  the  boundary 
layer.  Also  tested  was  the  technique  reported  by  Schmidt  and  Patankar 
(1988)  in  which  the  production  term  in  the  k  equation  is  set  to  zero  in 
laminar  regions. 

The  heat  transfer  distribution  predicted  by  setting  the  eddy 
viscosity  equal  to  zero  with  Tu=6%  is  shown  in  Figure  3.  The  results 
for  this  simulation  duplicate  the  convective  heat- load  distribution  from 
a  purely  laminar  simulation  (not  shown).  Results  from  setting  the 
production  of  A'  equal  to  zero  with  7w=6%,  shown  in  Figure  3,  indicate 
only  a  slight  increase  in  the  predicted  heat- load  levels  around  the 
airfoil  compared  to  the  prediction  with  0.  From  these  simulations 
it  is  concluded  that  neither  method  for  modeling  laminar  flow  with 
elevated  freestream  turbulence  accurately  predicts  the  wall  gradients  of 
the  quasi -laminar  boundary  layer. 

A  model  was  developed  based  on  studies  performed  with  the  data 
of  Arts  et  al.  (1 990).  Results  based  on  the  current  model  for  quasi  - 
laminar  boundary  layers,  labeled  as  4tQL  model”  (Praisner  and 
Rangwalla,  2001),  are  shown  in  Figure  4.  As  seen  in  this  figure,  the 
results  from  the  QL  model  are  more  accurate  in  quasi-laminar  regions 
than  the  results  shown  in  Figure  3.  The  current  method  is  based  on 
physical  reasoning  which  links  the  production  terms  in  the  k  and  o> 
equations  to  the  concept  of  self-sustaining  turbulence  in  turbulent 
boundary  layers.  The  analogy  is  drawn  that  in  laminar  regions  of  the 
boundary  layer,  where  disturbances  are  damped  by  the  mean  velocity 
profile,  the  production  of  both  k  and  to  should  be  negligible.  Implicit 
in  this  analogy  is  that  the  eddy  viscosity  in  a  quasi-laminar  boundary 
layer  is  independent  of  the  mean  strain.  Minimizing  the  production 
terms  in  the  k  and  &  equations,  in  contrast  to  setting  pf=  0,  allows  for 
the  convection  and  diffusion  of  freestream  turbulence  into  quasi- 
laminar  boundary  layers.  For  the  length  scales  (i.e.  dissipation  rates) 
and  turbulence  intensities  present  in  the  data  of  Arts  et  al.  (1990),  the 
QL  model  captures  convective  heat  loads  to  within  ±10%  in  quasi- 
laminar  regions  for  all  conditions  reported. 


Additional  testing  of  the  QL  model  was  performed  with  the 
cascade  data  of  Ames  (1994)  with  similar  accuracy  for  convective  heat 
loads  for  levels  of  Tu  up  to  1 2%  from  a  simulated  combustor  (Praisner 
et  al,  2003).  Additionally,  predictions  of  stagnation-point  heat  transfer 
levels  from  the  QL  model  are  within  ±10%  of  the  correlation  of  Van 
Fossen  et  al.  (1994)  for  a  wide  range  of  turbomachinery-specific 
turbulence  and  dissipation  levels. 

To  verify  that  the  quasi-laminar  boundary  layer  shape  is  also 
modeled  accurately  by  the  QL  model,  a  validation  study  was 
performed  with  the  favorable  pressure-gradient,  flat-plate  data  of  Blair 
(1981).  Figure  5  is  a  plot  of  nbn-dimensionalized  velocity  and 
temperature  profiles  in  the  quasi-laminar  flow  region  immediately 
upstream  of  transition  onset  with  ^=0.2x10'*  and  7w=3%  nominally, 
where  K  is  the  acceleration  parameter.  Also  shown  are  results  from 
CFD  predictions  performed  with  the  laminar  regions  simulated  by 
setting  ptj  -  0  (purely  laminar)  and  the  QL  model  Both  the 
momentum  and  thermal  boundary  layer  shapes  are  closely  matched 
with  the  QL  model  as  compared  to  the  method  of  setting  pT=  0.  Also, 
the  accuracy  of  the  predicted  convective  heat  transfer  levels  within  the 
quasi-laminar  region  were  as  accurately  captured  as  those  in  the 
benchmarking  against  the  data  of  Arts  et  al  (1990). 


Figure  5.  Non-dimensional  momentum  and  thermal  boundary 
layer  profiles  in  a  quasi-laminar  boundary  layer  just  prior  to 
transition.  Data  is  from  Blair  (1981). 


Figure  4.  Results  from  CFD  simulations  run  with  the  QL 
model  for  capturing  pre-transitional  quasi-laminar  boundary 
layers. 


For  the  momentum  boundary  layer  shown  in  Figure  5  the 
predicted  values  of  momentum-thickness  Reynolds  number,  Re&  are 
23%  and  6%  below  the  measured  value  for  the  pj=0  technique  and  QL 
model  respectively.  In  turn,  the  shape  factors  predicted  with  the  pj=0 
technique  and  QL  models  are  21%  and  3%  above  the  measured  level 
of  2.57,  respectively.  Also,  the  approximate  location  of  transition 
onset  for  the  case  with  Tu=6%  in  Figure  4  is  at  s=28mm  on  the  suction 
side.  The  QL-model  predicted  values  of  Re0  for  this  streamwise 
location  differ  by  1 5%  between  the  two  turbulence  levels  considered. 
These  results,  taken  together  with  the  flat-plate  results  for  Re&  indicate 
that  any  form  of  modeling  that  does  not  accurately  capture  the 
turbulence-enhanced  transfer  characteristics  of  quasi-laminar  boundaty 
layers,  will  provide  inaccurate  quantities  for  transition  onset 
prediction. 

Additionally,  if  one  assumes  the  Re $  value  for  the  pf*  0  technique 
is  essentially  an  exact  numerical  solution  for  a  condition  with  7w=G%, 
these  results  are  in  contrast  writh  the  assumption  made  by  Roach  and 
Brierley  (2000)  that  integral  quantities  of  pre-transitional  boundary 
layers  are  unaffected  by  freestream  turbulence.  The  aforementioned 
authors  supported  their  assumption  by  referencing  turbulent  boundary 
layer  data  with  elevated  Tu  levels.  Since  turbulence  levels  within 
turbulent  boundary  layers  are  of  the  order  10-15%  with  low  freestream 
Tu,  it  is  not  surprising  that  they  are  not  appreciably  affected  by 
elevated  levels  of  freestream  turbulence.  In  contrast,  the  low  levels  of 
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turbulence  present  in  laminar  boundary  layers  means  they  are 
susceptible  to  significant  momentum-  and  heat-transfer  enhancement 
with  the  addition  of  turbulence  convected  and  diffused  from  the 
freestream.  The  results  reported  here  are  consistent  with  the  apparent 
increases  in  measured  Re&  values  that  may  be  as  high  as  30%  in  pre- 
transitional  boundary  layers  subjected  to  elevated  levels  of  Tu 
measured  by  Yams  (2002). 

It  was  found  that  the  effects  of  roughness  on  pre-transitional 
boundary  layers  were  not  accurately  captured  with  the  combination  of 
the  QL  model  and  the  modified  wall  boundary  condition  for  m  from 
Wilcox  (1998)  for  rough  walls.  This  is  not  surprising  considering  that 
the  formulation  of  Wilcox's  modification  was  specific  to  turbulent 
boundary  layers.  Howev  er,  the  assumption  of  admissible  roughness  is 
typically  valid  for  all  but  the  high-pressure  turbine  in  large  commercial 
and  military  engines.  However,  in  the  case  of  small-scale  engines, 
such  as  those  used  for  unmanned  aircraft,  the  surface  roughness  levels 
from  cast  airfoils  may  play  an  important  role  in  the  transition  process* 
Resolution  of  pre-transitional  boundary  layers  subjected  to  significant 
levels  of  roughness  is  a  topic  that  requires  more  attention. 

Results  obtained  with  the  QL  model  for  various  values  of 
turbulent  Prandtl  number,  PrTt  show  little  variation  in  boundary  layer 
qualities.  Variation  of  PrT  from  a  molecular  level  of  0.7  to  a  turbulent 
level  of  0.9  resulted  in  variations  in  predicted  heat  flux  and  skin- 
friction  levels  of  less  than  5  percent  in  quasi -laminar  regions*  This 
may  be  related  to  the  fact  that  turbulent  mixing  in  quasi- laminar 
boundary  layers  is  significantly  lower  than  in  turbulent  boundary 
layers.  While  not  comprehensively  validated,  it  appears  that  there  may 
be  no  need  to  employ  different  turbulent  Prandtl  numbers  for  laminar 
and  turbulent  regions  of  the  flow  fie  Id  in  simulations  that  employ  the 
QL  model, 

MODELING  ATTACHED-FLOW  TRANSITION 

Having  demonstrated  an  ability  to  capture  the  development  of 
freestream  turbulence  level  and  dissipation  rate  as  well  as  pre- 
transitional  quasi- laminar  boundary  layer  details,  a  database  of 
experimental  test  cases  was  constructed  to  evaluate  methods  for 
transition  onset  prediction.  The  construction  involved  performing 
laminar  simulations  of  experimental  test  configurations  with 
appropriate  freestream  and  wall  boundary  conditions.  Inlet  levels  of 
turbulence  and  dissipation  rate  were  set  employing  the  aforementioned 
techniques  for  matching  the  decay  rate  of  Tu.  The  k~co- based  QL 
model  for  capturing  the  pre-transitional  boundary  layers  was  employed 
as  well.  For  each  experimental  case  considered,  the  measured  and 
predicted  static-pressure  distributions  were  in  close  agreement. 
Transition  onset  was  considered  to  occur  in  the  data  sets  where  wall- 
gradient  quantities  first  began  to  deviate  from  fully  laminar 
simulations.  Only  cascade  geometries  of  turbomachinery- specific 
airfoils  were  considered,  and  both  open-literature  and  in-house  data 
were  used  to  build  the  database  of  57  cases.  The  57  cases  themselves 
consisted  of  seven  different  geometries  tested  in  experiments  with 
various  boundary  conditions* 

The  database  does  not  necessarily  represent  accurate  physical 
quantities,  but  instead  represents  the  best  available  estimates  of  flow 
quantities  that  can  be  provided  by  two-equation  turbulence  model- 
based  RANS  simulations*  The  primary  use  of  the  database  was  to 
determine  if,  given  accurate  freestream  turbulence  and  pre-transitional 
boundary  layer  modeling,  accurate  transition-onset  predictions  could 
be  obtained  with  RANS  simulations  across  a  broad  range  of  flow 
conditions  and  geometries. 


A  Review  of  Models  for  the  Onset  of  Transition 

The  most  successful  physical  model  for  the  transition  process  is 
that  due  to  Emmons  (1951),  Schubauer  and  Klebanoff  (1955),  and 
Narasimha  (1957)*  In  this  model,  transition  is  considered  to  be  the 
result  of  the  random  formation  of  “spots"  of  turbulence  in  the 
boundary  layer  over  some  finite  region  in  the  streamwise  direction* 
These  turbulent  spots  grow  as  they  convect  downstream,  and  the 
intermitteney  (i.e.  the  fraction  of  time  the  flow  is  turbulent)  increases 
in  the  streamwise  direction  until  the  entire  surface  is  covered  with 
them.  At  that  point,  the  boundary  layer  is  considered  fully  turbulent. 

The  foregoing  wrould  suggest  that  an  assessment  of  the  effects  of 
various  flowfield  parameters  on  transition  must  be  an  evaluation  of  the 
way  in  which  they  affect  the  formation  and  subsequent  growth  of 
turbulent  spots  and/or  wake-induced  turbulent  strips.  Several  authors 
have  conducted  experimental  studies  in  that  regard  (e.g*  Clark  et  al, 
1994,  Gostelow  et  ah,  1994,  and  Halstead  et  al.,  1997),  and  others  (e*g. 
Chen  and  Thy  son,  1971,  Bey  and  Narasimha,  1990,  and  Solomon  et 
al*,  1996)  have  incorporated  information  on  turbulent-spot  formation 
and  kinematics  into  integral  methods  for  calculating  intermitteney  and, 
by  extension  through  a  “linear  combination  model"  after  Dhawan  and 
Narsimha  (1958),  skin  friction.  Still  others  (e.g*  Suzen  and  Huang, 
2000  and  Steelant  and  Dick,  2001)  have  made  use  of  the  “universal" 
intermitteney  distribution  to  derive  transport  equations  for 
intermitteney  that  are  solved  alongside  the  RANS  equations  or  used 
some  combination  of  the  correlations  of  Abu-Ghannam  and  Shawr 
(1980),  Drela  (1995),  and  Solomon  et  al*  (1996)  to  evaluate 
intermitteney  (e*g.  Gier  et  al.,  2000  Thermann  et  al*,  2001,  Roux  et  al*, 
2002,  and  Roberts  and  Yarns,  2003)  through  transition.  Subsequently, 
some  of  the  same  authors  modified  the  calculated  eddy  viscosity 
according  to  the  fraction  of  time  the  flow  is  predicted  to  be  turbulent 
through  the  transition  zone.  As  previously  stated,  pre-transitional 
boundary  layers  with  elevated  Tu  levels  may  not  be  adequately 
modeled  by  setting  jut  equal  to  zero  upstream  of  transition  onset* 

The  review  article  of  Mayle  (1991)  spurred  renewed  interest  in 
the  ideas  of  Emmons  (1951)  and  Narasimha  (1957)  for  the  prediction 
of  transition  in  turbomachines,  and  much  of  the  recent  work  described 
above  followed  recommendations  from  that  paper  closely.  In 
particular,  the  concept  of  the  “universal”  intermitteney  distribution  as 
described  in  detail  by  Narasimha  (1985)  has  been  used  to  build 
correlations  for  transition  onset  (Mayle,  1991),  turbulent-spot 
generation  rates  (Fraser  et  al.,  1994,  Gostelow  et  al.,  1994),  and 
transition  length  (Fraser  et  ah,  1994,  Gostelow  et  ah,  1994). 

At  the  root  of  many  correlations  used  to  predict  transition  onset  is 
the  F(t)  technique  of  Narasimha  (1985),  whereby  measured 
streamwise  variations  of  intermitteney,  %  are  plotted  in  the  form 

F(q  =  [-ln(l-r)  ]U1  (6) 

Narasimha  (1957)  proposed  that  all  turbulent  spots  are  formed 
randomly  in  time  and  span  wise  location  at  a  single  streamwise  station 
in  the  flow*  Under  that  hypothesis,  which  is  now  known  generally  as 
"concentrated  breakdown,"  Equation  6  becomes 

(7) 

for  constant-velocity  flow  along  a  flat  plate*  In  equation  7,  xt  is  the 
streamwise  location  where  all  spots  are  formed,  n  is  the  number  of  spots 
formed  at  that  position  per  unit  time  and  per  unit  spanwise  distance,  U„ 
is  the  freestream  velocity,  o  is  Emmons  (1951)  non-dimensional  spot- 
propagation  parameter,  and  x  is  a  location  on  the  fiat-plate  surface 
further  downstream  than  xt. 
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Dhawan  and  Narasimha  (1958)  also  showed  that  variations  of 
intermittency  through  the  transition  zone  possess  a  high  degree  of 
similarity  if  the  streamwise  distance  is  suitably  non-dimensionalized. 
They  demonstrated  that  the  transitional  intermittency- variations  from  a 
number  of  experiments  collapsed  on  a  single  curve  when  the  streamwise 
distance  was  represented  by 


i= 


x-xt 


(S> 


where  x15  and  x2$  are  the  streamwise  positions  where  the  intermittency  is 
0/75  and  0.25,  respectively.  The  same  authors  also  showed  that  the 
available  intermittency  data  were  well  represented  by  the  equation 

r  =  (9) 

Equation  9  has  become  known  as  the  " universal"  intermittency- 
distribution  (Narasimha,  1990  and  Narasimha,  1991).  Like  Equation  7} 
Equation  9  is  a  consequence  of  the  assumption  of  concentrated 
breakdown  as  applied  to  the  transition  model  of  Emmons  (1951)  in 
constant  velocity,  flat-plate  flow* 

Many  authors  have  shown  that  it  is  possible  to  linearize 
intermittency  distributions  using  Equations  6  and  7  and,  consequently,  to 
plot  the  data  in  the  universal  form  of  Equation  9  (e,g.  Owen,  1970),  This 
is  true  even  under  conditions  of  changing  pressure  gradient  and/or 
freestream  velocity  (Shanna  et  a L,  1982,  Gostelow  et  al.,  1994,  Fraser  et 
al.,  1994).  As  an  illustration  of  this  phenomenon,  representative  data 
from  Clark  (1993)  are  plotted  in  Figure  6.  Experimental  intermittency 
distributions  for  flat-plate  flows  under  both  favorable  and  adverse 
pressure  gradient  conditions  are  plotted,  and  both  agree  with  the 
universal  curve  very  well.  This  is  surprising  in  the  case  with  the 
favorable  pressure  gradient  since  the  local  velocity  varies  by  a  factor  of 
more  than  6.7  over  the  streamwise  distance  represented  in  the  figure. 
Note  that  both  the  spot  shape  parameter  and  die  freestream  velocity  itself 
are  considered  constant  and  brought  outside  an  integrand  in  the 
derivation  of  equation  7,  and  neither  of  these  assumptions  is  valid  in  the 
experiment  of  Clark  (1993).  Intermittency  variations  predicted  with  a 
time-marching  simulation  of  the  transition  zone  like  that  described  by 
Narasimha  (2003)  are  also  plotted  in  Figure  6.  These  predicted 
intermittency  distributions  are  taken  along  the  centerline  of  a  flat-plate 
flow  at  Mach  0.5.  Spot  propagation  parameters  were  as  measured  by 
Clark  et  al.  (1994),  and  three  different  distribution  functions  for  spot 
generation  were  considered.  The  universal  curve  fits  the  predicted 
intermittency  distributions  very  well,  not  just  when  concentrated 
breakdown  prevails,  but  also  for  both  a  bivariate  normal  (in  streamwise 
and  spanwise  directions)  spot  generation  function  and  a  point  source 
located  off  the  plate  centerline. 

Observations  like  those  presented  here  with  respect  to  the  universal 
curve  are  not  new  (See,  e.g.  Dhawan  and  Narasimha,  1958).  In  more 
recent  times  Mayle  (1991)  used  the  F{y)  technique  to  develop  a 
correlation  for  transition  onset,  where  onset  was  taken  to  be  the  x- 
intercept  of  the  F{y)  plot.  Also,  correlations  for  length  as  well  as 
turbulent- spot  generation  rates  under  a  variety  of  flow  conditions 
(Mayle,  1991,  Gostelow'  et  al.  (1994),  and  Fraser  et  al.,  1994)  have  been 
derived.  From  the  discussion  above,  it  does  not  follow  that  if  it  is 
possible  to  linearize  experimental  data  by  plotting  F(y ),  then  the 
transition  is  consistent  with  a  flat-plate  flow  undergoing  transition  via  a 
concentrated  breakdown  of  the  laminar  boundary  layer.  When  the 
current  database  for  attached-flow  transition  onset  is  compared  to  one 
such  correlation  from  Mayle  (1991),  as  in  Figure  7,  there  is  considerable 
scatter.  Testing  against  the  database  revealed  a  success  rate  of 
approximately  50%  in  predicting  onset  within  10%  of  the  measured 
location  in  terms  of  surface  distance.  Also,  transition  onset  typically 
occurs  at  lower  Reynolds  numbers  than  predicted  by  the  Mayle  (1991) 


relation.  Similar  results  to  those  in  Figure  7  are  reported  by  Simon  and 
Ashpis  (1996)  for  comparisons  made  between  data  and  die  Mayle  (1991) 
correlation.  In  light  of  these  findings,  it  might  be  better  to  develop 
correlations  for  transition  onset  and  turbulent-spot  production  rates,  for 
example,  from  direct  measurements  like  those  first  reported  by  Hofeldt 
(1996),  Such  direct  measurements  are  difficult,  however,  and  little  data 
is  available  at  present  to  create  such  correlations. 


Figure  6.  A  comparison  between  the  ‘‘universal”  curve  of 
Narasimha  (1985)  and  both  experimental  data  and  simulations 
from  Clark  (1993). 
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Figure  7,  A  comparison  between  the  current  database  for 
attached-flow  transition  onset  and  the  correlation  of  Mayle 
(1991). 

Both  Tani  (1969)  and  Reshotko  (1976)  review  a  number  of  models 
that  are  not  based  on  the  concept  of  universal  intermittency,  and  they 
both  point  out  that  one  of  the  first  was  due  to  Liepmann  (1945).  Many 
authors  refer  to  Liepmann  (1943)  as  the  source  of  this  idea,  but  the  1945 
publication  is  the  correct  one.  Liepmann  supposed  that  transition  onset 
occurs  when  the  local  Reynolds  stress  in  the  perturbed  laminar  boundary 
layer  equals  the  local  friction  velocity.  Liepmann" s  idea  has  been  recast 
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and  used  by  others.  For  example.  Van  Driest  and  Blumer  (1963) 
correlated  the  local  vorticity  Reynolds  number,  which  like  the  friction 
velocity  depends  on  the  normal  gradient  of  st  ream  wise  velocity  in  2D 
flow,  at  transition  onset  with  freestream  turbulence  and  pressure 
gradient.  Recently,  Mayle  and  Schulz  (1997)  referred  to  the  formulation 
of  Liepmann’ s  criterion  due  to  Shamia  et  al,  (1982)  as  appropriate  for 
transition  onset.  Sharma  et  al.  (1982)  argued  that  the  local  Reynolds 
stress  itself  depends  on  the  local  root  mean  square  of  stream  wise 
velocity  perturbations  and  stated  that  transition  onset  occurs  when 

u ' 

—  =  3.0  00) 

u 

based  on  experimental  data  related  to  a  turbine  airfoil  suction-side 
fiowfield.  In  equation  10,  u*  is  the  local  friction  velocity  and  u  '  is  the 
fluctuating  component  of  the  local  stream  wise  velocity.  Liepmann 
(1945)  argued  the  same  case  based  on  his  own  experimental  data. 
Follow  ing  Liepmann7  s  analysis  and  re-casting  his  correlation  in  the  form 
of  equation  10  results  in  a  constant  on  the  right  hand  side  equal  to  7.5 
rather  than  3.  In  addition,  Roach  and  Brierley  (2000)  report  that  the 
constant  in  equation  10  varies  from  3.0  to  7.4  for  the  ERCOFTAC  T3 
test  cases. 

If  one  assumes  that  ur  within  the  quasi-laminar  boundaty  layer  is 
modeled  to  a  reasonable  level  of  accuracy  with  the  QL  model,  it  might 
be  possible  to  predict  transition  onset  based  on  a  relation  like  that  in 
Equation  10.  In  Figure  8  the  quantity  u' !  u\  evaluated  at  transition 
onset,  is  plotted  for  all  cases  in  the  current  database.  For  comparison 
both  the  constants  of  Sharma  et  al.  (1982)  and  Liepmann  (1945)  are  also 
plotted.  Transition  onset  seems  to  occur  al  much  lower  levels  of  u  7  u 
than  those  indicated  by  Sharma  et  al.  (1982)  and  Liepmann  (1945)  over 
the  entire  database. 


Figure  8.  A  comparison  between  the  transition  onset  criteria  of 
Liepmann  (1945),  Sharma  et  al.  (1982),  and  Mayle  and  Schulz 
(1997)  and  the  current  attached-flow  database. 

The  poor  correlation  of  the  database  in  Figure  8  may  likely  be  a 
result  of  the  assumption  of  isotropic  turbulence  within  the  quasi-laminar 
boundary  layer  inherent  in  the  QL  model.  Also,  one  notes  that  the 
studies  of  Sharma  et  aL  (1982)  and  Liepmann  (1945)  pre-date  the  use  of 
conditional  sampling  techniques  for  boundary-layer  measurements  (See, 
e.g.,  Suder  et  al.,  1988  and  Kim  et  al.,  1994).  So,  it  could  be  that  the 
local  rms  of  velocity  fluctuations  in  the  boundary  layer  was  largely 
influenced  by  the  passage  of  turbulent  spots  over  the  fixed  hot  wires  in 
the  experiments  of  both  Liepmann  (1945)  and  Sharma  et  al.  (1982).  As 
such,  the  criteria  represented  by  equation  10  might  be  more  indicative  of 


turbulent-spot  detection  than  of  the  state  of  a  quasi-laminar  boundary 
layer  at  transition  onset.  The  foregoing  argument  is  supported  by  the 
results  of  a  simple  calculation.  Consider  a  boundary  layer  profile 
undergoing  transition  to  be  a  linear  combination  of  laminar  and  fully 
turbulent  profiles  weighted  on  the  intermittency,  as  suggested  first  by 
D  haw  an  and  Narasimha  (1958),  and  let  onset  occur  in  a  Blasius 
boundary  layer  on  a  flat  plate  at  a  local  streamwise  Reynolds  number  of 
250,000.  If  one  represents  the  turbulent  profiles  by  a  1/7 lh  povrer  law 
and  presumes  there  is  no  change  in  boundary  layer  thickness  at  transition 
onset  (strictly  speaking,  there  can  be  no  change  in  the  local  momentum 
thickness),  then  aty+  =  25,  the  local  rms  becomes  about  3  times  the  local 
friction  velocity  if  the  intermittency  is  10%.  This  is  very  much  in 
keeping  with  the  results  of  Sharma  et  al.  (1982),  In  any  regard,  the 
results  shown  in  Figure  8  suggest  that  even  with  the  current  modeling  of 
pre-transitional  boundary  layers,  accurate  RANS-based  transition 
predictions  can  not  be  obtained  with  a  model  of  the  type  in  Equation  10. 

A  correlation  for  transition  onset  that  was  developed 
independently  of  the  concept  of  intermittency  is  due  to  Abu-Ghannam 
and  Shaw  (1980).  Their  work  was  recast  by  Drela  (1995)  using  the 
boundary  layer  shape  factor  as  a  correlating  parameter  and  collectively 
these  relations  were  used  by  a  number  of  others  (e.g.  Fraser  et  al., 
1994,  Gier  et  ah,  2000,  and  Roux  et  al,  2002).  Abu-Ghannam  and 
Shaw  (1980)  measured  the  time- mean  near- wall  streamwise  velocity 
with  a  hot  wire,  and  they  determined  the  start  of  transition  to  occur 
where  that  parameter  began  to  deviate  from  the  reference  laminar 
level.  As  such,  their  method  is  consistent  with  demarcating  transition 
start  by  an  increase  in  the  local  skin  friction  over  the  laminar  level,  and 
that  is  in  essence  the  technique  used  to  determine  onset  for  the  current 
database.  The  percent  difference  between  the  momentum-thickness 
Reynolds  number  predicted  with  the  correlations  of  Abu-Ghannam  and 
Shaw  (1980)  and  the  actual  values  found  in  the  current  database  is 
plotted  in  Figure  9.  It  was  found  that  the  correlation  of  Abu-Ghannam 
and  Shaw  (1980)  produces  errors  in  the  momentum-thickness 
Reynolds  number  in  excess  of  20%  in  39  of  the  57  cases  in  the 
database.  Testing  against  the  database  revealed  a  success  rate  of 
approximately  68%  in  predicting  onset  within  10%  of  the  measured 
location  in  terms  of  surface  distance.  Testing  of  the  correlation  utilizing 
inlet  turbulence  levels  instead  of  local  levels,  which  is  consistent  with 
the  original  formulation  of  the  model,  did  not  improve  the  correlation 
of  the  database. 
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Figure  9.  A  comparison  of  the  efficacy  of  the  current 
transition-onset  model  for  attached-flow  with  that  of  Abu- 
Ghannam  and  Shaw  (1980). 
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A  Review  of  Models  for  Transition  Length 

Early  in  the  course  of  this  work,  three  separate  transition-length 
correlations  were  evaluated  in  conjunction  with  the  onset  model  of 
Abu-Ghannam  and  Shaw  (1980)  to  predict  the  development  of 
interm  ittency  through  transition  for  flow  fields  in  the  current  database. 
These  were  the  models  of  Gostelow  et  al.  (1994),  Fraser  et  al.  (1994), 
and  Solomon  et  ah  (1996).  However,  the  modeling  of  transition  length 
was  eventually  discarded,  and  abrupt  trips  at  onset  were  used 
exclusively  in  this  effort.  This  decision  was  made  for  a  number  of 
reasons.  First,  the  present  authors  thought  it  was  inadvisable  to  stack 
any  correlation  for  length  upon  another  correlation  for  transition  onset 
since  each  relation  must  have  its  own  attendant  uncertainty.  Further, 
the  models  of  Gostelow  et  ah  (1994)  and  Fraser  et  ah  (1994)  each 
assume  that  transition  length  is  a  function  simply  of  the  freestream 
turbulence  intensity  and  the  pressure  gradient  parameter  at  onset. 
Because  onset  on  a  turbine  airfoil  often  occurs  in  a  region  where  the 
pressure  gradient  changes  sign,  neither  length  model  was  found 
acceptable.  So,  the  model  of  Solomon  et  ah  (1996),  which  considers 
the  local  effects  of  pressure  gradient  on  turbulent  spot  growth  in  the 
determination  of  the  transition  length,  was  tested  Although  that 
model  was  expected  to  perform  better  on  turbine  airfoils  and  did  so, 
the  predictive  capability  of  the  relation  was  still  lacking.  Moreover,  it 
was  noted  that  all  three  length  models  depended  on  correlations  for 
turbulent-spot  production  rates  that  were  derived  from  experimental 
data  using  the  F(f)  technique.  On  account  of  the  discussion  presented 
above  relative  to  the  F(y)  technique,  none  of  these  models  were 
adopted.  Finally,  the  implementation  of  a  transition-length  model 
requires  that  a  streamline-based  application  of  the  modeling  be 
implemented.  That  is,  it  becomes  necessary  to  track  particles  from 
transition  onset  through  3D  flowfields  along  airfoil  surfaces.  So, 
without  the  use  of  an  intermittency-transport  equation  (Suzen  and 
Huang,  2000),  such  a  model  implementation  becomes  very 
cumbersome. 

Again,  based  on  the  difficulties  and  inaccuracies  inherent  in 
transition-length  modeling,  it  was  decided  only  to  model  transition 
onset  and  to  use  simple  instantaneous  (i.e.  abrupt)  trips.  Since  an 
abrupt  trip  implies  the  sudden  increase  in  turbulence  production  in  the 
context  of  this  work,  and  because  there  is  some  finite  distance  over 
which  the  boundary  layer  and  turbulence  model  respond  to  that,  a 
finite  transition  length  is  obtained  in  practice  (Praisner  et  al,  2003). 

A  New  Model  for  Attached-Flow  Transition 

Morkovin  (1978)  coined  the  term  bypass  transition,  and  has 
alerted  designers  to  the  implications  of  such  phenomena  many  times. 
In  his  NATO  AGARDograph  he  stated  that  designers  must  become 
“rather  sensitive  to  bypasses  when  transition  risks  imply  risks  to  the 
basic  mission  of  the  design,”  Such  a  situation  obtains  in  axial  flow 
turbomachinery  both  in  the  compressor  and  turbine  modules. 
Morkovin  also  pointed  out  that  “it  is  easy  to  criticize;  it  is  another 
matter  to  offer  a  constructive  suggestion  to  the  designer.”  With  these 
ideas  in  mind  and  with  the  current  capability  to  predict  local  flow 
features  at  the  onset  of  transition  for  experiments  pertinent  to 
turbomachinery  design,  the  database  was  analyzed  to  build  a  new 
model  for  attached-flowr  transition  onset. 

It  was  postulated  that  the  boundary-layer  momentum  thickness  at 
the  onset  of  transition  depends  upon  a  number  of  local  flow 
variables  at  the  edge  of  the  boundary  layer  and  the  local  wall 
temperature  (TV).  These  flow'  variables  were  taken  to  include  the 
density  (p),  the  dynamic  viscosity  (//),  the  sonic  speed  (tf),  the  driving 
temperature  for  heat  transfer  (Tg\  the  thermal  conductivity  of  the  fluid 
(ft),  the  specific  heat  at  constant  pressure  (c^),  the  magnitude  of  the 
flow  velocity  ({/„),  the  streamwise  velocity  gradient  (dC/^Ms),  the 


root- mean-square  of  streamwise  velocity  fluctuations  (w  %  and  their 
length  scale  (>£). 

Applying  the  Buckingham-Pi  Theorem  via  the  “step-by-step” 
technique  of  Massey  (1986),  it  is  possible  to  represent  the  physical 
process  as  a  function  of  7  non-dimensional  parameters.  The  relation  is 

Ree=  f[M,Tg  ITW,  K,Tu,  Pr,  A/0  )  00 

The  non-dimensional  parameters  in  equation  1 1  are  the  momentum- 
thickness  Reynolds  number  at  transition  onset  ( Re# ),  the  Mach  number 
(M),  the  gas-to-wall  temperature  ratio  (Tg  /TV),  the  acceleration 
parameter  (K)>  the  local  turbulence  intensity  (Tu\  the  molecular 
Prandtl  number  (TV),  and  the  ratio  of  turbulent  length  scale  to 
momentum  thickness  (A  /&).  Note  that  the  physical  dimensions  used  to 
define  the  indicial  matrix  are  mass  (m),  length  (£),  time  (f), 
temperature  (I),  and  heat  (H).  Because  heat  is  invoked  as  a 
dimension,  it  is  assumed  that  there  is  no  substantial  conversion  of 
mechanical  energy  into  thermal  energy.  So,  the  results  are  restricted  to 
flows  of  gases  with  moderately  supersonic  Mach  numbers,  but  this  is 
not  a  limitation  for  the  gas-turbine  situation. 

The  non-dimensional  parameters  in  Equation  1 1  were  tabulated 
for  each  case  in  the  database,  and  the  results  were  correlated  to 
determine  the  combination  of  non-dimensional  parameters  that 
provided  the  best  collapse  of  the  data.  Tire  best  correlation  for 
transition  onset  was  found  when 

Tuii  “2> 

where  A  and  B  are  constants  equal  to  8.52  and  -0.956,  respectively. 
The  entire  database  is  plotted  in  Figure  10.  Equation  12  implies  that 
the  local  momentum  thickness  at  transition  onset  is  a  function  of  only 
one  non-dimensional  parameter  that  is  a  product  of  two  of  the  basic 
non-dimensional  parameters  derived  directly  from  the  “step-by-step” 
method  of  Massey  (1986). 


Tu(0)/2, 

Figure  1 0.  The  current  model  for  the  onset  of  attached-flow 
transition  compared  to  the  database. 

An  ability  to  determine  the  edge  of  the  boundary  layer  in  a  robust 
manner  is  an  important  aspect  of  transition  modeling  techniques  based 
on  integral  and  edge  quantities.  The  edge  of  the  boundary  layer  was 
taken  to  be  the  distance  from  the  wall  at  which  the  local  vorticity 
dropped  to  1  %  of  the  maximum  value  in  the  O-grid  at  that  streamwise 
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grid  location.  This  method,  reported  by  Michelassi  et  al.  (1998), 
provides  a  robust  and  accurate  edge  detection  technique  for  laminar 
boundary  layers.  Variations  in  the  percentage  of  the  maximum 
vorticity  used  to  define  the  boundary  layer  edge  between  0,8%  and 
1.8%  resulted  in  variations  in  the  coefficient  of  correlation  of  the  least- 
squares  fit  for  Equation  12  of  between  0.95  and  0.96.  The  boundary 
layer  grid  densities  outlined  in  the  Computational  Methods  section 
were  necessary  for  the  implementation  of  the  vorticity-based  edge 
detection  method. 

As  seen  in  table  1,  all  the  non-dimensional  parameters,  save  the 
molecular  Prandtl  number,  varied  markedly  at  transition  onset  over 
this  database.  For  example,  transition  onset  occurred  under  both 
favorable  and  adverse  pressure  gradients.  The  strength  of  the  pressure 
gradients  covers  the  range  from  above  relaminarization  (AT  >  3x10**) 
on  the  favorable  side  through  Thwaite's  separation  criterion  (A^  Re#2  < 
-0.09,  See  White,  1991)  under  decelerating  conditions.  Both  adiabatic 
flows  and  those  with  velocity  profiles  that  were  affected  by  the 
temperature  dependence  of  viscosity  were  part  of  the  database.  Also, 
the  range  of  Mach  numbers  covered  incompressible  through  transonic 
flows.  The  local  turbulence  intensity  varied  by  two  orders  of 
magnitude,  and  while  the  maximum  value  may  at  first  consideration 
seem  low,  one  must  remember  that  the  in  lei  turbulence  intensity  was 
much  higher.  Further,  the  ratio  of  local  turbulence  length  scale  to 
momentum  thickness  at  onset  varied  by  an  order  of  magnitude.  Again, 
the  molecular  Prandtl  number  was  essentially  constant  for  all  points  in 
the  database. 

Table  1 ;  Variation  of  non-dimensional  parameters  at  transition  onset. 


Variable 

Range 

Re# 

73  -  856 

Kx  IQ6 

-1.9  -  4.8 

KRee2 

-0.15  -  0.057 

M 

0.05  -  1.24 

Tg/Tw 

L0  -  1.41 

TUotmet  (%) 

0.1 1  -  5.09 

A/e 

4.26  -  66.2 

Pr 

0,71  -  0.71 

Physical  Significance  of  the  Current  Attached-Flow  Model 

Equation  12  was  incorporated  into  the  airfoil  design  system 
directly  with  the  factor  and  power  in  the  equation  set  equal  to  8.52  and 
-0.956,  respectively.  However,  it  is  possible  to  recast  the  relation  in 
such  a  way  that  gives  insight  into  its  physical  significance.  One  notes 
that  the  constant  B  in  the  relation  is  very  close  to  -I.  If  that  value  is 
accepted,  then,  after  some  algebraic  manipulation.  Equation  12 
becomes 


where  A}  is  another  constant  that  may  be  evaluated  directly  as  the 
mean  of  values  occurring  at  transition  onset  in  the  database.  Here,  A  t 
was  found  to  be  7.0  ±  1.1.  Note  that  Equation  13  implies  that 
transition  onset  occurs  wrhen  the  ratio  of  a  boundary-layer  diffusion 
time  {p  0  2/p)  to  a  timescale  associated  with  the  large-eddy  turbulent 
fluctuations  ( le=/Ju ' )  becomes  a  critical  value. 

It  is  instructive  to  consider  the  implications  of  Equation  1 3  for  a 
Blast  us  boundary  layer  undergoing  transition.  The  usual  form  of  the 
laminar  diffusion  timescale,  is  p  5  2/p  (See  e.g.,  Schlichting,  1 979, 
and  Hofeldt  et  al.,  1998)  where  5  is  the  thickness  where  the  local 
velocity  becomes  99%  of  the  freestream  value  and  p  is  the  fluid 
density.  The  ratio  of  99%  velocity  thickness  to  momentum  thickness 


in  a  Blasius  boundary  layer  is  5:0.664.  Substituting  into  Equation  13, 
taking  A  j  equal  to  7.0,  and  rearranging  gives 

—  0-25  td  (14) 

The  time  for  a  Blasius  boundary  layer  to  grow  to  a  thickness  Sis  1  /25th 
the  laminar  diffusion  timescale  (Schlichting,  1979).  So,  transition 
onset  would  occur  for  the  Blasius  boundary  layer  when  the  local  eddy 
timescale  approached  a  level  associated  with  a  growth  of  ^6  boundary 
layer  thicknesses.  Schlichting  (1979)  gives  the  smallest  unstable 
wavelength  of  a  Blasius  boundary  layer  as  =6£  So,  Equation  14 
implies  that  the  onset  of  bypass  transition  occurs  when  the  local  eddy 
timescale  reaches  a  timescale  associated  with  the  wavelength  of  a 
Tol  lmi  en-Schlichting  (TS)  wave. 

Although  one  typically  does  not  associate  TS  activity  with  bypass 
transition,  it  has  been  noted  by  Herbert  (1988)  that  the  ultimate 
breakdown  associated  with  the  appearance  of  “spikes”  in  hot-wire 
records  and  the  first  formation  of  turbulent  spots  in  natural  transition 
occurs  over  a  length  scale  of  about  1  TS  wavelength.  Also,  Walker 
and  Gostelow  (1990)  have  measured  frequencies  consistent  with  TS 
activity  in  boundary  layers  undergoing  bypass  transition  in  adverse 
pressure  gradients.  Additionally,  Mack  (1975)  noted  that  often  TS 
frequencies  persist  in  boundary  layers  undergoing  natural  transition 
beyond  the  range  of  applicability  of  small  disturbance  theory.  Volino 
(2002)  also  reported  that  for  separated-flow  transition,  TS  frequencies 
were  detected  for  both  lowr  and  high  levels  of  7n.  So  perhaps  the 
results  of  this  exercise  are  not  so  surprising.  Again,  the  real 
implication  of  this  discussion  is  that  when  the  ratio  of  the  turbulent- 
eddy  time  scale  to  the  laminar  diffusion  time  reaches  a  critical  value, 
bypass  transition  occurs.  Further,  the  critical  value  of  this  ratio  is 
nearly  constant  over  a  range  of  flow  conditions  consistent  with  gas- 
turbine  engines, 

MODELING  LAMINAR  SEPARATION  AND  TURBULENT 
REATTACHMENT 

“Of  all  the  transition  modes,  there  is  none  more  crucial  to 
compressor  and  low-pressure  turbine  design  and  none  more  neglected 
than  separated-flow  transition”  (Mayle,  1991).  As  attempts  are  made 
to  reduce  airfoil  counts,  and  hence  component  cost  and  weight,  airfoil 
loadings  need  to  increase.  Highly  loaded  airfoils  are  more  prone  to 
experiencing  laminar  separations  (Figure  11)  and  stall.  In  the  stalled 
condition  profile  losses  can  increase  as  much  as  500%  over  the  case 
where  the  separation  reattaches  to  the  airfoil.  Laminar  separations 
occur  in  the  leading-edge,  suction- side  regions  of  compressors  and  the 
aft  suction-side  regions  of  Low-Pressure  Turbine  (LPT)  airfoils. 
Airfoil  stall  causes  compressor  surge  and  poor  performance  of  LPT s  at 
cruise  conditions. 

If  the  laminar  shear  layer  formed  by  a  separation  transitions  to  a 
turbulent  state  close  enough  (i.e.  in  a  stream- wise  sense)  to  the 
separation  location,  it  typically  reattaches  to  the  airfoil  surface  as  a 
result  of  turbulent  mixing  that  entrains  high-momentum  fluid  into  the 
near-wall  region.  This  scenario  is  schematically  depicted  in  Figure 
11a.  However,  if  transition  of  the  shear  layer  occurs  sufficiently  far 
downstream  of  separation,  the  layer  typically  does  not  reattach, 
resulting  in  a  stalled  condition  as  shown  in  Figure  lib.  It  is  therefore 
critical  that  a  design  system  be  capable  of  predicting  the  existence  of  a 
laminar  separation  and  whether  or  not  the  laminar  separation  will 
reattach. 
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Figure  11.  Schematic  representation  of  suction-side,  laminar- 
separation  characteristics  showing  both  reattached  (a)  and 
stalled  (b)  conditions. 

As  a  first  step  in  developing  a  transition  modeling  capability  for 
separated  flow,  a  proof-of-concept  study  was  executed  with  the  goal  of 
determining  if  a  separated-flow  transition  model  could  be  effectively 
implemented  in  a  RANS  code.  For  this  study,  CFD  simulations  of 
cascade  experiments  were  run  with  imposed  point- like  trips  set 
according  to  experimental  hot-film  data.  Figure  12  is  a  plot  of 
separation  and  reattaehment  locations  as  a  function  of  exit  Reynolds 
number  for  a  cascade  airfoil  with  7w=5%  (data  if  from  Butler  et  al., 
1988).  The  total  suction-side  surface  length  of  the  airfoil  for  this  case 
was  37.1  mm.  Experimental  and  computational  results  are  shown. 


Figure  12.  Measured  and  predicted  separation  and 
reattachment  locations.  Transition  was  specified  in  the 
simulations  based  on  data. 

It  can  be  seen  in  Figure  12  that  the  separation  and  reattachment 
locations  from  simulations  utilizing  the  QL  model  for  the  pre¬ 
transit  ional  boundary  layers  are  in  close  agreement  with  the  data.  The 
deviation  from  the  data  at  the  lowest  Reynolds  number  was  a  result  of 
stall  occurring  in  the  simulation  and  subsequent  unsteadiness  of  the 
separation  location.  For  a  Reynolds  number  of  1.5xl06  a  simulation 
was  also  run  with  the  turbulent  viscosity  set  to  zero  in  laminar  regions 
and  the  results  are  also  plotted  in  Figure  12.  For  this  simulation  the 
separation  location  occurs  upstream  of  the  measured  location.  This  is 
consistent  with  the  quasi-lam  mar  boundary  layer  shape,  and  hence 


near- wall  momentum,  not  being  accurately  modeled  by  setting  /*t=0. 
Similar  comparisons  were  made  with  additional  cascade  data  at  a 
variety  of  turbulence  levels  and  the  results  were  consistent  with  those 
in  Figure  12.  The  calculated  loss  levels  from  the  simulations  shown  in 
Figure  12  were  in  vety  close  agreement  with  the  measured  values. 
The  conclusion  of  this  study  was  that  the  RANS  code  with  abrupt  trips 
and  the  QL  model  for  pre-transitional/separation  boundary  layers 
would  provide  an  accurate  framework  for  separated-flow  transition 
modeling. 

So,  in  addition  to  the  attached-flow  transition  studies,  a  CFD- 
supplemented  database  of  laminar- separations  with  turbulent- 
reattachment  was  constructed  based  on  47  in-house  and  open -literature 
experimental  cascade  data  sets.  The  test  cases  for  the  separated-flow 
database  include  laminar  separations  with  turbulent  reattaehment  on 
both  compressor-  and  turbine-specific  geometries.  Some  of  the  data 
sets  included  in  the  database  are  from  airfoils  on  the  verge  of  stall. 
Like  the  attached-flow  database,  the  separated-flow  database  covers  a 
significant  range  of  turbomachinery-specific  flow  parameters. 

A  Review  of  Models  for  the  Onset  of  Sep  a  rated -Flow 
Transition 

One  common  method  for  predicting  the  transition  location  of  a 
near- wall  bounded  shear  layer  involves  the  use  of  a  correlation  that 
was  developed  for  attached- flow  transition  such  as  either  the  Mayle 
(1991)  or  Abu-Ghannam  and  Shaw  (1980)  model  These  models  rely 
on  a  non-dimensional  boundary  layer  thickness  (Re#  or  shape  factor) 
for  the  prediction  of  transition  onset  and  were  developed  with  data  for 
attached- flow  transition.  As  depicted  in  Figure  11,  in  the  separated 
region,  the  only  attached  boundary  layer  that  exists  is  the  reverse-flow 
boundary  layer  within  the  separation  bubble.  The  primary  issue  with 
employing  an  attached-flow  transition  model  in  a  region  where  the 
flow  is  separated  is  that  the  physical  significance  of  Re#  fundamentally 
changes  once  the  boundary  layer  separates. 


IF  =  Tu(C  IX)02 


Figure  13.  A  comparison  between  the  separated-flow  transition 
model  of  Roberts  (1980)  and  the  separated-flow  transition 
database. 

A  rather  unique  model  for  the  prediction  of  separated-flow 
transition  was  reported  by  Roberts  (1980),  in  which  the  distance  from 
the  separation  location  to  the  shear  layer  transition  location  (L  in 
Figure  1 1)  is  related  to  freestream  turbulence  quantities.  This  model  is 
unique  in  that  it  considers  both  turbulence  intensity  and  length  scale  in 
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predicting  separated-flow  transition  onset.  Figure  13  is  a  comparison 
between  predictions  performed  with  the  model  of  Roberts  (1980)  and 
the  separated- flow  database.  Plotted  in  this  Figure  is  Reynolds 
number  based  on  L  versus  the  local  ‘Turbulence  Factor”  which  is 
defined  as  TF=Tu  (C  fX  )0^,  where  C  is  the  airfoil  chord.  The  highest 
level  of  turbulence  factor  considered  in  the  development  of  the  Roberts 
(1980)  model  was  approximately  0.06.  This  model  was  developed  for 
external  flows  with  low  ffee-stream  turbulence  levels  (<0.2%).  The 
only  database  cases  that  are  well  correlated  by  this  model  are  the  cases 
in  which  7w<Q,6%  (Figure  13).  The  Roberts  model  (1980)  does  not 
correlate  database  cases  with  Tu>0,6%  well  enough  for 
implementation  in  a  design  system.  Additionally,  alteration  of  the 
model  by  leaving  out  the  length  scale  in  the  calculation  of  77%  as 
suggested  as  a  possible  modification  by  Roberts  (1980),  did  not 
improve  the  correlation  of  the  database. 

Other  separated-flow  transition  models  have  been  reported  by 
Walker  (1989),  Mayle  (1991),  and  Batman  and  Wang  (1999).  These 
models  relate  separation  length  to  the  conditions  of  the  laminar 
boundary  layer  at  the  separation  location.  RAN S -based  simulations 
employing  models  of  this  type  have  proven  to  be  at  least  trend  accurate 
in  the  prediction  of  separated  transition  (Volino,  2002,  Houtennans  et 
al.,  2003).  The  separated-flow  transition  database  was  employed  to 
test  the  models  of  Mayle  (1991)  for  “long”  and  “short”  bubbles  and  a 
comparison  of  the  predicted  and  measured  bubble  lengths  is  shown  in 
Figure  14.  It  can  be  seen  in  this  figure  that  neither  the  long-  nor  the 
short-bubble  model  provides  sufficient  accuracy  for  design  purposes. 
Similar  results  were  obtained  for  comparisons  between  the  models  due 
to  Walker  (1989)  and  Batman  and  Wang  (1999)  and  the  separated- 
flow  database.  These  results  are  supported  by  Volino  (2002)  who 
reported  that  the  correlations  of  Batman  and  Wang  (1999),  Mayle 
(1991),  and  Davis  et  al.  (1985)  give  “rough”  estimates  for  the 
transition  behavior  of  his  experimental  test  cases. 


‘Separation! 

Figure  14.  A  comparison  between  the  separated-flow  transition 
models  of  Mayle  (1991)  and  the  separated-flow  transition 
database. 

A  New  Model  for  Separated-Flow  Transition 

In  the  context  of  RANS  simulations  the  authors  believe  that  the 
multi -moded  nature  of  transition  is  best  captured  by  employing  two 
separate  models  for  attached  and  separated  transition.  So,  following 
the  body  of  material  concerning  models  for  separated-flow  transition 
summarized  above,  a  new  model  has  been  developed  based  on  the 


separated-flow  database.  The  same  dimensional-analysis  technique 
used  for  the  attached- flow  model  was  employed  in  the  development  of 
this  model.  The  best  correlation  of  the  database  was  obtained  when 
the  length  of  the  bubble  was  related  to  the  state  of  the  boundary  layer 
at  separation.  The  form  of  the  current  model  is: 

~  =  CRc0.sepD  (15) 

where  C  and  D  are  constants  equal  to  173.0  and  -1.227,  respectively,  L 
is  the  distance  between  separation  and  transition  onset  and  Siep  is  the 
surface  distance  from  the  stagnation  point  to  the  separation  location. 

The  separated-flow  database  is  shown  in  Figure  15  in  the  context 
of  the  current  transition  model.  The  reasonably  good  correlation  of  the 
database  with  the  recasting  of  Walker's  original  idea  supports  the 
assertion  that  for  viscously  dominated,  near- wall  bounded  separations, 
the  bubble  size  scales  on  the  state  of  the  boundary  layer  at  separation. 
While  Equation  15  does  not  explicitly  contain  turbulence  quantities, 
they  are  still  important  in  the  determination  of  Re0  at  the  separation 
location  if  a  model  such  as  the  current  QL  model,  is  employed  to 
capture  the  effects  of  freestream  turbulence  on  the  pre-separation 
quasi -laminar  boundary  layer.  Additionally,  it  should  be  noted  that 
equation  15  was  implemented  in  the  design  system  in  a  conservative 
fashion  by  setting  C  at  a  level  20%  higher  than  the  least- squares -fit 
value  plotted  on  Figure  1 5  as  “current  model.”  This  conservatively 
defined  correlation  is  also  plotted  in  Figure  1 5  as  “shifted  model.” 


Figure  15.  The  current  model  and  database  for  separated-flow 
transition.  The  model  with  a  conservative  shift  is  also  shown. 


In  his  experimental  assessment  of  the  Pack  B  airfoil,  Volino 
(2002)  reported  that  boundary  layer  reattachment  occurs  essentially 
simultaneously  with  transition  onset  for  separated  flow  transition. 
This  supports  the  assertion  made  by  Lou  and  Hourmouziadis  (2000) 
that  the  transition  length  is  very  short  because  the  lack  of  wall 
damping  in  the  shear  layer.  Also,  Walker  et  al.  (1988)  reported  that 
abrupt  transition  is  a  relatively  (compared  to  transition-length  models 
for  separated  flow)  realistic  model  of  the  transition  region  in  laminar 
separation  bubbles.  So,  as  in  the  case  of  attached-flow  transition, 
abrupt  transition  is  assumed  for  the  implementation  of  the  current 
separated-flow  model  in  the  RANS  solver.  Validation  of  the  many 
assumptions  involved  in  the  development  of  both  models  as  well  as 
their  implementation  is  reported  in  Part  II  of  this  work. 
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Physical  Significance  of  the  Current  Separated-Flow  Model 

While  the  concept  of  “long”  and  “short”  bubbles  has  been 
discussed  in  literature  relating  to  separated- flow  transition  for  many 
decades,  the  present  results  suggest  that  the  concept  may  not  be 
important  in  turbomachinery  configurations.  The  present  database  is 
in  conflict  with  reports  of  Hatman  and  Wang  (1999)  and  Houtermans 
et  al  (2003)  where  the  existence  of  two  bubble  regimes  seems  evident 
from  data.  In  the  current  database,  no  “long”  bubbles,  that 
significantly  alter  airfoil  pressure  distributions,  were  found  that  also 
reattached  to  the  airfoil  In  other  words,  separations  that  do  not 
reattach,  and  hence  alter  the  turning  characteristics  and  loadings  of 
airfoils,  are  not  referred  to  as  bubbles  here  as  they  are  not  closed. 
Large  reattaching  bubbles  that  significantly  alter  pressure  distributions 
are  possible  on  flat-plate  configurations  such  as  that  employed  by 
Hatman  and  Wang  (1999)  if  reattachment  occurs  downstream  of  the 
simulated  trailing  edge.  Therefore,  their  correlation  may  still  be  useful 
in  predicting  where  the  laminar  shear  layer  transitions  downstream  of 
an  airfoil  trailing  edge.  The  present  database  shows  a  continuous 
distribution  of  separation  behavior  up  to  the  point  of  stall  This  is 
consistent  with  the  results  of  Houtermans  et  al.  (2003)  for  their  “short” 
bubble  regime.  The  data  regarded  as  reflecting  “long”  bubble 
behavior  by  Houtermans  et  al  (2003)  is  here  interpreted  as  arising 
from  stalled  conditions.  So,  the  authors  believe  that  the  “short”  bubble 
regime  primarily  describes  separation  behavior  with  reattachment.  In 
this  work,  if  separated- flow  transition  onset  is  predicted  to  occur 
downstream  of  the  airfoil  trailing  edge,  then  the  airfoil  is  run  fully 
laminar. 

CONCLUSIONS 

The  ability  to  model  accurately  the  development  of  free  stream 
disturbances  has  been  shown  to  be  an  important  aspect  of  any 
transition  modeling  capability.  Evidence  has  been  presented  that  the 
effects  of  freestream  turbulence  on  pre-transitional  boundary  layers 
should  be  accurately  modeled  before  any  form  of  empiricism  is 
applied  for  the  prediction  of  transition.  The  k-m  model,  along  with  the 
current  method  for  modeling  quasi -laminar  boundary  layers,  was 
employed  to  supplement  104  cascade  data  sets  from  literature  and  in- 
house  studies  (i.e.  Pratt  &  Whitney  proprietary  data)  to  build  databases 
for  attached-  and  separated-flow  transition.  Existing  models  for  both 
transition  mechanisms  were  assessed  with  the  databases  and  deemed 
not  to  have  sufficient  accuracy  for  design  purposes.  Consequently, 
dimensional  analysis  was  employed  as  a  guide  for  the  extraction  of 
pertinent  flow  variables  from  flowfield  predictions  of  the  experiments, 
and  new  models  were  developed. 

Two  correlations  were  developed,  one  for  the  onset  of  attached- 
flow'  transition  and  the  other  for  the  length  of  a  separation  bubble  prior 
to  turbulent-reattachment  The  current  models  are  based  on  local 
flowfield  parameters,  and  they  appear  to  have  greater  efficacy  than  a 
number  of  extant  correlations.  In  particular,  the  model  for  attached- 
flow  transition  appears  to  have  a  physical  basis  with  respect  to  the 
fundamental  mechanism  of  bypass  transition  in  compressible  flow. 
That  is,  it  was  found  that  the  onset  of  transition  occurs  when  the  ratio 
of  a  boundaty-layer  diffusion  time  to  a  timescale  associated  with  the 
local,  energy -bearing  turbulent  fluctuations  at  the  edge  of  the  shear 
layer  reaches  a  critical  value.  Further,  it  was  found  that  the  critical 
value  of  the  ratio  was  nearly  constant  over  a  wide  range  of  flowfield 
conditions  consistent  with  turbomachinery  airfoils.  By  contrast,  no 
such  underlying  physical  basis  was  apparent  from  considerations  of 
the  separated-flow  model:  it  appeared  to  be  more  of  a  straight 
correlation  of  variables. 

The  models  have  been  implemented  as  point- wise  trips  in  a  3D 
RANS  solver  that  forms  part  of  a  turbomachinery  design-system.  It 


should  be  noted  that  both  the  attached-  and  separated-flow  models  are 
based  solely  on  two-dimensional  data  and  applications  of  them  in 
three-dimensional  flow  fields  may  elucidate  deficiencies.  In  addition, 
no  modeling  has  been  implemented  to  account  for  the  effects  of 
roughness  on  pre-transitional/separation  boundary  layers.  Part  II  of 
this  paper  focuses  on  the  validation  of  the  current  models  for  use  in  an 
airfoil  design  system 
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